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Abstract 

Well-known and newly developed renormalization schemes for tan (3 are an- 
alyzed in view of three desirable properties: gauge independence, process in- 
dependence, and numerical stability in perturbation theory. Arguments are pro- 
vided that no scheme can meet all three requirements, and as an illustration, a 
"No-Go-Theorem" for the renormalization of tan (3 is established. Nevertheless, 
two particularly attractive schemes emerge. A discussion about which scheme 
might be the best compromise in practice is given. 

PACS numbers: 12.60.Jv, ll.lO.Gh, 11.15.-q 

1 Introduction 

The quantity tan /3 is one of the main input parameters of the minimal supersymmetric 
standard model (MSSM). At the tree level it is defined as the ratio of the two vacuum 
expectation values f 1^2 of the MSSM Higgs doublets, 

ian 13=—. (1) 

Vl 

Owing to its central appearance in spontaneous symmetry breaking, tan /? plays a cru- 
cial role in almost all sectors of the MSSM and has significant impact on most MSSM 
observables. However, the vacuum expectation values and tan /3 are not directly mea- 
surable quantities. In spite of its importance, tan [5 is an auxiliary variable. The virtue 
of tan [5 is that it can be used as an easy-to-handle input parameter in terms of which 
all different observables can be expressed. In contrast with other parameters, such as 
the electron charge or the particle masses, there is no obvious and unique way to relate 
tan (3 to an observable. 

The actual definition of tan /3, its physical meaning and its relation to observables 
is given by the choice of a renormalization scheme. This choice determines the nu- 
merical value of tan (3 as well as its formal properties, such as gauge dependence and 
renormalization-scale dependence. 
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In the literature, several renormalization schemes for tan f3 have been proposed 
and used [1-3]. Each of these schemes has specific advantages and disadvantages. In 
this paper, we give three criteria that are desirable for a renormalization scheme for 
tan [3: 

• Gauge independence. If the relation between tan f3 and observables is gauge 
independent, the numerical value of tan (3 is also gauge independent and a more 
physical interpretation of tan (3 is possible. 

• Process independence. If tan (3 is defined by a relation to a physical process, 
a non-universality and flavor dependence can be introduced that violates the 
intuition implied by eq. O that tan /? is a quantity of the MSSM Higgs sector. 
Process-dependent schemes have further technical drawbacks like the necessity 
to calculate more complicated vertex functions in the determination of 5 tan /? 
and the possible appearance of infrared divergent QED or QCD corrections. 

• Numerical stability. The numerical properties of the renormalization constant 
6 tan P should not spoil the validity of the perturbative expansion. Generally 
speaking, the finite contribution to 6 tan /3 and the renormalization-scale depen- 
dence of tan (3 should not be too large. 

We analyze the properties of the known renormalization schemes in view of these crite- 
ria and devise new schemes in order to satisfy them. Since it is the most intricate point, 
we will be concerned with questions of gauge dependence and gauge independence for 
the largest part of the paper. 

Our results are negative. We find that the known process-independent schemes 
are gauge dependent and that the new gauge-independent schemes lead to numerical 
instabilities. Arguments are provided that schemes satisfying all three criteria do not 
exist. However, as a result of our analysis two schemes emerge that are particularly 
attractive and useful compromises. 

The outline of the present paper is as follows. After briefly introducing the MSSM 
Higgs sector in Sec.|2l it is shown in Sec.|3]that the DR scheme as well as the schemes 
proposed in [1,2] are gauge dependent already at the one-loop level. In Sec.|4j three 
gauge-independent and process-independent schemes are developed. The numerical 
instability induced by these schemes is exhibited in Sec.|5] The discussion of process- 
independent schemes is completed in Sec. |6l by demonstrating that a large class of 
process- and gauge-independent schemes leads to numerical instabilities. Finally, in 
Sec.0the problems caused by process-dependent schemes are discussed and one use- 
ful scheme is presented. Section [8l contains our conclusions and a discussion of two 
attractive renormalization schemes for tan (3. 

2 The MSSM Higgs sector 

The MSSM contains two Higgs doublets 

rr _( ^l + j-,i<Pl-^Pl)\ rr_( \ 
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whose electrically neutral components are shifted by f i 2 in order to account for the 
finite vacuum expectation values. The parameters of the MSSM Higgs sector are 



9,9','ml,rhl,rnl, 



(3) 



where g and g' are the SU(2) and U(l) gauge couplings, 771^2 = /^^ + "^12 ^^^^ 
the Higgsino parameter /x, and mf 2 3 are soft-breaking parameters. We assume CP 
conservation, so all parameters are real. Expressed in terms of these parameters, the 
Higgs boson potential reads 



V = mf|i?ir + m^\H2r + ml{H{H^ - HfH^ + h.c.) 



(4) 



In order to make the physical content more transparent, the parameters of the gauge 
and Higgs sector ^ and the shifts vi^2 are reparametrized in terms of 



e, Mz, Mw, Ma, tan /3, ti,t2 
with {g'^ = 9^ + g'"^, = vl+ V2) 



99' ,.2 fv'^ 
e = -^, Mz = iv.^ 

9 ^ 



M, 



2 9'^v'^ V2 

~ -rr^ tan/3 = — , 
2 vi 



t2 = V2 



\f2v\ \plv2 
mfvi + mlv2 + -^viivf - f|) 

mlv2 + mlvi - —V2{vl - vl) 



m|(tan (3 + cot /?), 



(5) 

(6a) 
(6b) 

(6c) 

(6d) 



For vanishing tadpole parameters ti = t2 = 0, M\ is the mass of the pseudoscalar A^, 
and the tree-level Higgs mass eigenstates are given by the fields h, G^,A^, , 
with 
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C/3 
S/3 

C/3 
S/3 



-S/3 
Cf3 

-S/3 
C/3 



and 



tan 2a = tan 2(3 



M\ + Ml 



Ml 



H 

h 



-Tr/2 < a < 0. 



(7) 
(8) 
(9) 

(10) 



Here and in the following we use the abbreviations tp = tan (3, s/3 
cos P, etc. 



sin f3, C/3 
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The tree-level Higgs-boson masses read 

M^Jl^ =Ml + M^, (11) 
=\(^Ml + M'z± ^(M2+M|)2-4M|M2cos2(2/3)) . (12) 

An important equation is given by 

jy^2(0)j^2 (0) 

''"''^^^^ " M2(M2W+Mf(0)-Ml)- ^^^^ 

It relates tan /3 to the physical Higgs-boson masses and will be useful for later pur- 
poses. 

The counterterms to the MSSM Higgs sector are generated by a multiplicative 
renormalization transformation of the parameters and fields. For the purpose of the 
present paper, the explicit form of the field renormalization is not important, but the 
parameter renormalization is given by 

ml ^ ml + 6ml, fi,2 ^ ^^1,2 + 5^1,2- (14) 

Through the relations ^ also renormalization constants 6e, 6Mz,w,A, ^tp and 5ti^2 
are defined, in particular 61/3 = tan/3(^ ~ l^^)- There are two important relations 
between 61^ and other renormalization constants (for ti = t2 = G) that will be of use 
later: 

— + — = V2{6ml + 6ml) + V26ml{tp + ^) + ^ml{l - \)6tp, (15) 
6Ml = -6ml{^ + ^) - - ^)6t, + 4-|i- + c}^^. (16) 



3 Gauge dependence of some schemes for tan /3 
3.1 Definition of the schemes 

The actual definition of the parameters their physical meaning and their relations 
to experimental quantities is given by the choice of a renormalization scheme. In 
gauge theories a gauge fixing is necessary for quantization. Since the gauge fixing 
is unphysical, the relations between observable quantities do not depend on it, but 
the relations between observables and the parameters ^ can be gauge dependent. 
Accordingly, in such a case the values for the parameters extracted from experiment 
are gauge dependent. 

For e, Mz^w,A^ on-shell renormalization scheme (see [4] for the case of the 
MSSM) provides a gauge-independent definition [5]: e is related to the effective charge 
in the Thomson limit, Mz^w,A are the masses of the Z, W, bosons (defined as the 
real parts of the poles of the respective propagators). In addition, setting the renormal- 
ized tadpoles to zero, 

r^, = r^, = 0, (17) 
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is in agreement with the requirement of gauge independence. Here and henceforth, 
r denotes the generating functional of the renormaUzed IPI vertex functions and 
^ipi - {Pij • • •) denotes a vertex function with incoming fields (pi . . . and incoming mo- 
menta pi 

In contrast, for tan (3 no such gauge-independent standard-definition is available. 

In this section we will show explicitly the gauge dependence of several well-known 
schemes for tan/?. The first of these is the DR scheme, which is defined by the 
condition 



DR : Stp = pure divergence, (18) 

where "pure divergence" denotes a term of the order A = — Je + log 47r in di- 
mensional reduction (or in dimensional regularization when suitable supersymmetry- 
restoring counterterms have been added [6]) whose prefactor is such that all renor- 
malized quantities are finite. Further common renormalization schemes for tan /3 are 
the ones introduced by Dabelstein [1] and by Chankowski et al. [2] (DCPR). In these 
schemes one writes 6vi^2 = t'i,2(^-^i,2/2 — Svi^2 with the field renormahzation con- 
stants 2 of the Higgs doublets and requires the conditions 

— = —, ReS^02 Ml = 0, (19) 

where the renormalized A'^Z two-point function is decomposed as ^AOzt^(.~P:P) = 
ip^t^Aozip'^)- At the one-loop level these conditions lead to 

DCPR : Stfs = ^^^Re^AOziMl) (20) 

with the unrenormalized A^Z two-point function S^o z- A similar prescription requir- 
ing vanishing H'^W~ mixing instead of vanishing A^Z mixing has been used in [7]. 



3.2 Extended Slavnov-Taylor identity as a tool 

In the presented schemes tan /? is apparently not directly related to any observable. 
It is known that the DR scheme leads to a gauge dependent tan (3 at the two-loop 
level [8]. In the following, we will show that in fact both schemes lead to a gauge 
dependence of tan /3 already at the one-loop level. 

The tool we use to determine the gauge dependence is an extended Slavnov-Taylor 
identity introduced in [9, 10]: 

S{T) = S{T) + xd^T = 0. (21) 

Here ^ denotes an arbitrary gauge parameter in the gauge-fixing term and x is a 
fermionic variable acting as the BRS transformation of ^. S{T) is the usual Slavnov- 
Taylor operator (see Appendix). If not stated otherwise we leave the form of the gauge- 
fixing term open. We only assume that the gauge-fixing term is coupled to auxiliary B 
fields, 

L&,= J2 B^'^J^'' + ^-\B''\\ (22) 
V=A,Z,W+,W- 
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This simplifies the symmetry identities, but for loop calculations, eq. (l22t is equivalent 
to Lfix = — ^ X] l-^^P' which is obtained by eliminating the B fields using their 
equations of motion. 

In contrast to the usual Slavnov-Taylor identity S{T) = 0, the extended identity 
(ETl need not be satisfied. If, however, S{r) = holds, then physical quantities 
— expressed as functions of the parameters of the lowest-order Lagrangian — are 
gauge independent. Accordingly, S{r) = implies that the values for the parameters 
extracted from experiment are gauge independent. If suitable x-dependent terms are 
added to the action, the extended Slavnov-Taylor identity holds at the tree level. In 
order to satisfy S{T) = also at higher orders, the renormalization conditions must 
be chosen so as to not contradict S{T) = 0. 

We will make use of the extended Slavnov-Taylor identity in particular in two 
ways. On the one hand, in our one-loop calculations we employ regularization by 
dimensional reduction [1 1] assuming all symmetries are preserved, so S^(r(^^'''^s) = 0, 
where r'''^^ denotes the unrenormalized vertex functional and the index denotes 
the loop order. The identity ^^(r^^)'''''^) = g entails an easy way to calculate the ^- 
derivatives 9^r^^.^.'.!°^ of regularized vertex functions. 

On the other hand, presuming S{T^'^^) = at some loop order, there is a very sim- 
ple necessary condition for ^(r) = at this order, yielding a simple check whether 
a given set of renormalization conditions is compatible with S{T) = 0. As shown 
in [10, 12], the counterterms in gauge theories can be generally classified into gauge- 
dependent total BRS variations and gauge-independent non-BRS variations. Indepen- 
dent of the presence of spontaneous symmetry breaking, the non-BRS variations corre- 
spond to the symmetric parameters of the theory, whereas the vacuum expectation val- 
ues as well as field renormalization constants are related to total BRS variations. This 
result can be easily transferred to the case of the MSSM Higgs sector (see also [4] 
for the classification of the MSSM counterterms). It means that in order to satisfy 
S{T) = 0, 

d/rSg^'^ = d^Smj^ = d^Smj = (23) 

has to hold. In contrast, for 6vi^2 and field renormalization constants a ^ dependence 
is compatible with ^(r) = 0. Equation (l23l also constitutes a sufficient condition for 
S{T) = 0, since if it is satisfied, x-dependent counterterms can always be added in 
such a way that S{r) = holds at the considered loop order. 

Thus, eq. d23l provides a simple way to check whether a given renormalization 
scheme is compatible with S(r) = 0. In the following calculations we will use the 
on-shell conditions for e, Mz,w,A together with the tadpole conditions dTTl . Since 
these conditions are compatible with S{r) = 0, eq. d23l provides in particular a check 
for the renormalization condition for tan P and hence for the gauge independence of 
tan/3. 

3.3 Calculating the gauge dependence 

Let us now show the gauge dependence of tan /? at the one-loop level in the DR and 
DCPR schemes (CHl, beginning with the 7JR scheme. If S{T) = is to hold, the 
gauge independence of the counterterms (l23l imphes, together with eq. (fTSl . a certain 
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gauge dependence of the renormalization constant 5tg 



fin. 



— H 



V2ml ( 1 - ^ 



(24) 

Vi V2 



We only consider the purely finite parts here since the divergent contributions are re- 
stricted by 5(r(^)''^'^s) = and hence are in agreement with S{T) = 0. In the DR 
scheme the finite part of the renormalization constant 6tp itself is zero and therefore 
gauge independent, 

Se^t/-^ = 0, (25) 

and the question is whether this is compatible with the gauge dependence prescribed 
by eq. d2^ . 

The tadpole counterterms in (l24l are determined by F,^^ j = 



Sh,2 = r(^)f ^. (26) 

Using the extended Slavnov-Taylor identity at the regularized level, 5^(r(^)'''®s) = q, 
yields for the l.h.s. of (l24t : 

fSh St2\ sr^ p(l),reg v;^ J_p{0) 

^ ^ ^ ' j=l,2 i=l,2 

where the Y^^ denote the sources of the BRS transformations of (pi used in the Slavnov- 
Taylor identity. The ^ dependence of F^^^''^'^^ and the results for '^^^y^'^^ depend on 
the specific choice of the gauge fixing. In the gauge we obtain (for the relevant 
Feynman rules see the Appendix): 

^{i),reg _ / cos/3 \ 1 f gMz p ^,^2 ^,.^2^ 

+ ^i?o(0,eM^,eM2.)), (28) 
where Bq denotes the usual two-point function. Hence, in the context of the R^ gauge, 

a,(^ + ^)" = 0, (29) 

and the DR scheme is compatible with eqs. d25l and hence with S{T) = at the 
one-loop level. This confirms the result of [8], where one-loop gauge independence, 
but two-loop gauge dependence of tan /3 was found in the R^ gauge. 

However, in more general gauges (l29t is not true. As a simple example we consider 
an infinitesimal deviation from the R^ gauge by introducing a second, infinitesimal 
gauge parameter that enters the gauge-fixing function for the Z boson as follows: 

J^^ = df.Z'' + Mzi^G^ + C^^°A^). (30) 
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Figure 1: The one-loop diagrams to the vertex functions V^)/^'^^ and F 

X^<t>i X Ys, 



We can study the dependence on this second gauge parameter in the same way as the 
dependence on ^. We introduce the variable x^^° consider the Slavnov-Taylor 



identity S{T) = S{T) + x 



ZA" 



d^ZAO F 



0. As in the general case, the on-shell 



conditions for e, Mz,w,A and the tadpole conditions (fTTl are in agreement with S{T) = 
0.' Using this identity we obtain an identity prescribing the C^^" dependence of Jt/?^"^, 
analogously to eq. 



d^ZAO 



(6t,_ 




V ^1 


V2 ) 



V2m 







where, analogously to eq. d27t . the l.h.s. is given by 



(31) 



d^ZAO 



6ti 6t2 

Vi V2 



ml 



1 

^2 



_l_ „, -p(l),reg 
+ ^zaO 



Y 



02 



(32) 



However, the one-loop results of the vertex functions involving ' are (evaluated at 



c 



0) 



— sin (3 
cos (3 



1 gMz 



167r2 2 



Bo{0,^MlMl) 



Hence, the l.h.s. of eq. dSTT i is non-vanishing: 



6ti 6t^, 



d^zAO ( \ 



Vl 



V2 



fin 



(33) 



(34) 



Thus, in the generalized gauge, the (^^^° independence of tan /? and S{T) 

d^ZAO Stp''^ / 0. 



imply 
(35) 



'clearly, <30> induces a gauge-dependent mixing of A^' with and B^. One might wonder 

whether this mixing has an effect on the on-shell condition for Ma and whether it is still in agreement 
with the requirement of gauge independence. As shown in [4], as long as the auxiliary B fields are not 
eliminated from the gauge fixing <22t . this on-shell condition can be written as 



det 



(j? = xi) = 



with ReA^^ = M^. Moreover, with uneliminated B fields there is no contribution of the gauge fix- 
ing to rgOGO, Tgoyio, r^o^o, and the evaluation of the renormalization conditions at the tree level is 
unchanged compared to the 7?^ case. At higher orders, as shown in [4], this condition for Ma indeed 
corresponds to the (gauge-independent) pole of the A" propagator, even though the ^4" field mixes with 
B^ at the tree level — or with G" and Z'^ after the elimination of B^. 
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Figure 2: The one-loop diagrams to the vertex functions F^^oy^ and F^^oy^^. 

This is incompatible with the DiZ-condition, which implies d^zAoStf^^'^ = 0. Con- 
versely, if the DR scheme is used, the relation between tan /3 and observable quanti- 
ties, i.e. the value tan /3°^p extracted from experiment is gauge dependent: 

a^;,^otan/3|^/0. (36) 

The different result for eq. d33b compared to eq. d28l can be understood by considering 
the corresponding Feynman diagrams in Fig.^Cthe Feynman rules can be read off from 
the Lagrangian given in the Appendix). In the gauge fixing only the Golds tone 
bosons G^'^ appear, and correspondingly only Goldstone bosons contribute as internal 
scalar lines in the diagrams for F^y^ . In the generalized gauge, however, we study the 
dependence on the part of the gauge fixing, and correspondingly only appears 
as an internal scalar line in F^^^Oy^ . The couplings of G^'^ and A^ to {^^) are 

proportional to (g°^^) and ("qs"/?^), respectively, which explains the results d33l and 
(EH. 

In a second step it is quite easy to see the gauge dependence of tan /? in the DCPR 
schemes defined by eq. d20t . We know that the DR scheme is in agreement with 
5(F) = in the i?^ gauge at the one-loop level. So we consider the difference 

d^{6tf3BCPR - StpTm) ^ ImF^')^^g'fi'^(Ml). (37) 
Using S(r^^^) = 0, in the R^ gauge the r.h.s. can be expressed as 



i=l,2 ^=Z.G^ 



fin 

+ 0, (38) 



which is non-zero as can be easily seen by inspection of the corresponding one-loop 
diagrams (see Fig.|2li. For example, there is an M/i-dependent contribution to ^xA'^Yqc 
that cannot be cancelled in eq. ( I38t . 

Hence, in the DCPR schemes tan fi is gauge dependent even at the one-loop level 
in the i?^ gauge: 

5^tan/?g^c%^0- (39) 

4 Three gauge-independent schemes for tan /3 

In the following we present three gauge-independent renormalization schemes for 
tan /3. These schemes share the additional property that tan /? is defined via quan- 
tities in the Higgs sector and without reference to a specific physical process — in that 
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sense they are closely related to the tree-level definition tan P = The use of these 
schemes is on the one hand to demonstrate the existence of such gauge-independent 
schemes and on the other hand to illustrate three ways to devise gauge-independent 
renormalization schemes. 

The first scheme is defined by the requirement that the parameter m| and its coun- 
terterm is gauge independent as dictated by eq. d23l . which is most easily realized by 
DR renormalization: 

m-s scheme : 5m^ = pure divergence. (40) 

We refer to this scheme as the "7713 scheme". Together with the on-shell conditions for 
e, Mz,w,A, ti,2 this condition fixes the seventh parameter in tan (3. Since at the 
regularized level S'(r(^)''^°s) = holds, the divergent part of dm"^ is gauge indepen- 
dent, hence 

d^6ml = (41) 



in the ma scheme, and thus this scheme is compatible with eq. and S{T) = 0. 
Hence, it defines tan /? in a gauge-independent way. Using the relation (l6bt between 



Mj,, tan f3 and m| we can derive a result for the finite part of 6t 



A-f fill jtj. fin 

5Mlfi° = sin2/?^ + cos2/?^-mi(l-cot2/3)5t/°. (42) 
v2wi v2w2 

In this way Stfj is expressed via (5M^ and 6t 1^2 by the self energy and the tadpole 
contributions. The superscript denotes the purely finite part of the renormalization 
constants. The result for the purely divergent part of Stf^ is the same as in all other 
schemes. 

A second gauge-independent scheme can be read off from eq. (l24l . which is a 
necessary consequence of S{T) = 0. The most straightforward way to satisfy eq. (l24l 
is to require 

Tadpole scheme : = — + — ) • (43) 

V2m|(l - T^) V ^1 V2 J 

This scheme will be referred to as the "Tadpole scheme". It can also be obtained in a 
second, instructive way by temporarily replacing the tadpole conditions T^^ ^ = Qhy 
the conditions 

(5ti,2 = 0. (44) 

In such a scheme, where the tadpoles are not renormalized, the renormalization con- 
stants 5vi^2 are gauge independent if S(r) = holds (since in this case 5vi^2 can be 
entirely expressed in terms of the gauge-independent constants in eq. dl^). Hence 
also 5tp is gauge independent, and in combination with eq. d44l i the Di?-condition 

6tf/''' = (45) 

is in agreement with S{T) = 0. 
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The connected Green functions and the physical content of the theory do not 
change by varying 6vi^2 and accordingly 5ti,2, Stp, etc. while keeping 6g^'\ 5mf 2' 
(5m| fixed. Equation (fTSl yields a relation between 6ti 2, Stp and these fixed renor- 
malization constants: 

5*/"^ = ^ ' M f- + + ^(^<2,Sml) (46) 

V2m|(l - 1 ) V V2 J 

Therefore, for physical quantities the conditions F^^ ^ = for the tadpoles and eq. 
for tan /3 are equivalent to eqs. (1?^ and ( R3t . 

By construction, the Tadpole scheme (l43l combined with the on-shell conditions 
for e, VK,A5 *i,2 is in agreement with the extended Slavnov-Taylor identity SiV) = 
0. Hence it provides another gauge-independent definition of tan (3. 

In the third scheme, tan (3 is defined as a combination of physical Higgs-boson 
masses in agreement with the lowest-order result (fT3l : 

HiggsMass scheme : cos^(2/3) = , , , ^ , ttoy, (47) 

^ ^ M\{M% + Ml- M\) 

where Mfj ^ denote the physical masses of H, h (defined as the real parts of the poles 
of the respective propagators). This scheme is denoted as "HiggsMass scheme". Since 
the physical masses are gauge-independent quantities, this scheme provides manifestly 
a gauge-independent definition of tan /3. 

The definition (l47l is problematic, because due to higher-order corrections the 
r.h.s. can get larger than unity. In spite of this possibility, evaluating the requirement 
(l47l at the one-loop level yields a condition for 5tp that always has a solution: 

„ _!_ Sfc/i ^HH _ ^hh + '^HH 

where Tjhh,HH denote the renormalized hh and HH self energies, evaluated on-shell 
(j)^ = jj, respectively). This equation is linear in 6tf} and can therefore always be 
solved for dtp as a function of the unrenormalized self energies, 5M\ and Sti^2- 

However, the possibility of the r.h.s. in eq. ( HTl getting larger than unity already 
signals that this HiggsMass scheme could cause numerical problems. In the next 
section, the numerical properties of all three gauge-independent schemes will be dis- 
cussed. 



5 Numerical instability of the gauge-independent schemes 

Gauge dependence and gauge independence are important conceptual characteristics 
of renormalization schemes. But in practice it is also mandatory that a renormalization 
scheme does not lead to numerical instabilities in the quantum corrections to physical 
processes. The gauge-dependent schemes presented in section |3] have been success- 
fully used in many practical loop calculations. Only in the calculation of the neutral 
Higgs-boson masses [13] slight numerical instabilities caused by the scheme of [1] 



II 







tan P = 
large /i 


3 

no mixing 




tan (3 = 
large fi 


50 

no mixing 


DR 
Tadpole 


-0.06 
0.81 
4.50 


-0.06 
-0.46 
-0.21 


-0.06 
-0.04 
1.24 


-0.17 
285.29 
370.73 


-0.17 
127.11 
140.11 


-0.17 
4.92 
34.53 



Table 1: 9 tan /3/(?log jl in the DR- m^, and Tadpole scheme for various parameter 
scenarios. We have chosen Ma = SOOGeV, and the remaining parameter values are 
chosen according to [14]. 



were reported that are avoided using the DR scheme. In this section we will demon- 
strate that the three gauge-independent schemes presented in the foregoing section lead 
to much more severe numerical problems at the one-loop level. 

A straightforward way to analyze the numerical behavior of the one-loop cor- 
rections is to derive the renormalization-scale dependence of tan/? in the various 
schemes. In the DCPR schemes d20t and in the HiggsMass scheme Wfl . tan /3 is 
scale independent because it is defined on-shell, while in the DR scheme, the m-s 
scheme and the Tadpole scheme the scale-dependence of tan /3 is obtained by simply 
equating 

(tan/3 + 5tp{fl))^^p^ = (tan/3(/i) + <5t^(A))other scheme ' (49) 

where jl is the renormalization scale used in dimensional reduction. For the scale- 
dependence in the DR scheme we obtain a simple analytical formula, which is well- 
known (see e.g. [8]): 

— ^(tan/3(/Z))^ = tan(3-^{3hl-3h^t+hl), (50) 
alog/x Idtt^ ^ 

where ht^b,T are the Yukawa couplings of the top, bottom, and r, respectively, and 
where the contributions of the first two generations are neglected. The analytical re- 
sults for the scale dependence in the other schemes are more complicated, so we restrict 
ourselves to a numerical analysis. Tabled contains the numerical results for the scale 
dependence for three typical sets of MSSM-parameters taken from [14]. 

Apparently, the scale dependence of tan f3 in the DR scheme is quite small, indi- 
cating a small uncertainty in one-loop corrections due to 6tf}. On the other hand, the 
scale dependence of tan (3 in the ms and Tadpole schemes can get immensely large. 
In practice, for instance changing the scale fl from rrit to Ma can cause unaccept- 
ably large changes in the numerical values of one-loop corrections to observables that 
depend on tan /3. 

The fact that tan f3 is not scale-dependent in the HiggsMass scheme does not imply 
that this scheme does not lead to numerical instabilities. In table|2lthe numerical values 
of the purely finite part dtf^^'^ are shown in the schemes where it is non-vanishing. 
Typically, the finite contribution to a renormalization constant should be suppressed 
compared to the respective tree-level parameter. Table|2lshows that this is the case only 
for the DCPR schemes (and of course for the DR scheme), whereas in the three gauge- 
independent schemes, in particular in the HiggsMass scheme, 6tfj^^ can exceed tan /? 
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tan/? = 


3 




tan /3 = 50 






''''h 


large 


no mixing 


''''h 


large fi 


no mixing 


DCPR 


-0.10 


-0.06 


-0.08 


3.56 


UAl 


0.46 


m3 


0.56 


-0.08 


-0.04 


490.45 


-61. U 


-4.85 


Tadpole 


2.64 


-0.46 


0.33 


624.70 


-16.A6 


0.92 


HiggsMass 


-2.44 


-1.83 


-1.33 


-426.54 


-1995.93 


-314.50 



Table 2: (5t/3 in the DCPR, ms, Tadpole, and HiggsMass scheme for various param- 
eter scenarios. We have chosen Jl = Ma = 500GeV, and the remaining parameter 
values are chosen according to [14]. 



by far. The consequences are very large one-loop corrections to quantities depending 
on 5tp, signalizing the breakdown of the validity of the perturbative expansion. 

As an example for the influence on the calculation of observables we consider the 
one-loop results for the mass of the light CP-even Higgs boson, M^, according to the 
strict one-loop formula 

Ml = M2(0)-ReS^,(M2(0)), (51) 

where 'Ehh is the renormalized one-loop /i-self energy. The numerical results are pre- 
sented in table |3l for the case of small tan where the Higgs-boson mass is most 
sensitive to tan /3. 

We observe that for the parameters used here the DCPR schemes and the DR 
scheme agree well numerically, whereas the discrepancy between these schemes and 
the gauge-independent ones can be very large. These discrepancies cannot be inter- 
preted as theoretical errors in the prediction of the Higgs-boson mass but they are 
consequences of the invalidity of the perturbative expansion in the m^, Tadpole and 
HiggsMass scheme. 

In contrast, perturbation theory is trustworthy in the DCPR schemes and the DR 
scheme. A more detailed comparison between the scheme of [1] and the DR scheme 
taking into account the leading two-loop effects has shown that the DR scheme has 
generally a better numerical behavior in certain regions of the parameter space [13]. 

6 Impossibility of defining tan /3 in the Higgs sector without 
introducing gauge dependence or numerical instability 

In the previous section severe numerical problems were found in all three gauge- 
independent schemes. In this section it is shown that the conflict between gauge in- 
dependence and numerical stability is unavoidable if tan f3 is defined via quantities of 
the Higgs sector. Since the numerical instability is so strong that the schemes can- 
not be used in practice, this result can be interpreted as a "No-Go"-like theorem for a 
gauge-independent renormalization of tan f3 in the Higgs sector. 

More precisely, we study the class of renormalization schemes where 5t/3^'^ at the 
one-loop level is given by a linear combination of finite parts of the on-shell Higgs self 
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large fi 


no mixing 


Tree level 


72.51 


72.51 


72.51 


DCPR 


134.44 


97.40 


112.23 


DR, fl = MA 


135.03 


97.93 


112.83 


DR, fi = nit 


134.63 


97.38 


112.35 


HiggsMass 


119.58 


81.48 


102.87 


ms, ft = Ma 


138.34 


97.28 


112.56 


ms, fi = mt 


143.23 


93.23 


112.25 


Tadpole, jl = Ma 


149.94 


94.08 


115.18 


Tadpole, p, = mt 


173.45 


92.21 


123.97 



Table 3: The light Higgs mass Mh in the five renormalization schemes and for var- 
ious parameter scenarios. We have chosen tan (3 = 3 and Ma = 500GeV, and the 
remaining parameter values are chosen according to [14]. 



energies, their momentum derivatives S', and tadpoles: 

= linear combination of (^S^o^o(M|), I;^ogo(M|), S^o^(M|), T^HniM^), 

(52) 

The coefficients in this linear combination should be functions of the parameters of the 
Higgs potential e, tan /3, Mz,w,A- The choice of d52t is motivated by the intuition that 
tan /3 is a quantity of the MSSM Higgs sector. All schemes considered in sections 13 
|4]belong to this class. 

As shown in Appendix O the most general gauge-independent schemes of the 
class ( 1521 are given by 

dt^fi'^ = 6tfs^^ + (uaKa + aHKH + ahKh) (53) 

where aAM,h are coefficients and the quantities KA,H,h are defined as the following 
combinations: 

Ka = ^AOAO - E ^AOAO^. (rS)'.' (54) 

and analogous for KH,h- These combinations are gauge independent as can be seen 
from the identities (displayed here for the case of Ka) 

a X+ — T(l)-regp(0) o _ -p(l),reg-p(0) 

derived from S'(r(i)'''^s) = 0. Besides, Ka. H,h are nothing but the gauge-independent 
(5M^ ^ mass counterterms in the scheme with 5ti 2 = (compare discussion of the 
Tadpole scheme). 
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The numerical instability originates from several terms contributing to 6tp and the 
scale dependence c?iog tan /3. In particular, in the general scheme d53t . 9iog;itan/3 
contains terms of the order 

^ Ml 
Ml' Ml ' Ml' ^ ^ 

which can get large independently of each other. ^ It is possible to choose the coef- 
ficients aA.H,h in eq- 'EU such that these three large terms are exactly cancelled in 
(?iog p, tan l3. In this way an optimal, gauge-independent scheme with minimal scale 
dependence is obtained. Of course, this optimal scheme is nothing but the HiggsMass 
scheme d47l . where tan (3 is defined on-shell and therefore //-independent. As seen in 
section |5j the HiggsMass scheme leads to numerically not acceptable, large contribu- 
tions to 5tp itself. 

Hence, the only scheme where the contributions (l56l are absent is the numerically 
inacceptable HiggsMass scheme, and all other schemes of the form d53l involve the 
large contributions (l56l to ^log^tan/? and are for this reason unsuited. None of the 
gauge-independent schemes defined in eq. (l53l can be used in practice. 



7 Process-dependent schemes 

If tan (5 is defined via quantities of the Higgs sector there is an unavoidable conflict 
between gauge independence and numerical stability. In order to circumvent these 
problems one could try to define tan j3 outside of the Higgs sector by relating it to a 
specific physical process. This method has been adopted in [3], where it was suggested 
to use the decay —>■ r^Ur- The one-loop corrected decay width to this process 
reads 

amlMl^tl r 



2— + + 2^ - - 

e TJlr tp M^ sw 



(57) 



where Fu+^^ is the form factor describing the vertex and external wave-function cor- 
rections to the amplitude —>■ t^Vj- including the H^-W^ and H^-G^ mixing 
self energies. 

By requiring that the radiatively corrected decay width retains the same form as the 
lowest-order formula, eq. d57l can be understood as a definition of the renormalization 
constant 5tj^. As a consequence of the relation to a physical observable, this definition 
of tan /? is manifestly gauge independent. 

However, this scheme also has several drawbacks. At first, for the computation of 
5tp it is necessary to compute loop corrections to the three-particle vertex in Ffj+rv, 
which can be difficult beyond the one-loop level. Furthermore, it is conceptually dis- 
advantageous to define tan /? in a specific process, since in this way it becomes a 
non-universal, flavor-dependent quantity. 

^In many of the cases in tab. Q|2||3| these terms are actually subdominant due to the large value of 
Ma- Nevertheless, in general they are very significant, especially if the ratios in eq. <56> get large. 
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Finally, the decay vertex —>■ r+z/,- also receives QED corrections, which nec- 
essarily include contributions with real photon emission in order to cancel infrared 
divergences. It is not possible to separate the QED corrections from the rest of the 
electroweak corrections since they are not individually UV finite. For practical cal- 
culations, however, it is unacceptable to include real bremsstrahlung corrections into 
the definition of a counterterm, since this procedure should depend on experimental 
phase-space cuts and is therefore technically very involved. 

While the first two drawbacks hold for any process-dependent scheme, the problem 
posed by the QED corrections can be avoided by the choice of another process. One 
possibility is given by the decay A'^ —>■ t^t^. Its one-loop decay width reads 



T ' T 



1 + F^o, 



+ 2— + 2^ + 2^ - - 2^ 

e nir tp sw 



(58) 



The QED corrections to this decay width consist of the photon loop contributions to 

Fj^Ott and drriT, 



T-iQED , j-iweak 



(59) 



These QED corrections F'^f^ + 2 *^"^""^ form a UV-finite subset of the full elec- 
troweak one-loop corrections. The reason for this difference to the decay T~^i/.r 
is that the latter process relies substantially on the SU(2) symmetry and 7, Z and W 
loops have to be summed to yield a UV-finite result. In contrast, the QED corrections 
to — > r'^r" can be thought of as being generated by an effective theory containing 
essentially the A^tt vertex and the photon and are therefore naturally finite. 

Owing to the finite QED corrections, a possible definition of tan (3 is given by 
requiring that the pure weak corrections in eq. (l58l cancel and that the exact decay 
width is given by the tree-level result plus QED corrections: 



l + FX^:!^ + 2^ + 2^ 
e rrir 



weak 



T[A° 



T ' T 



+ 2^ 



5M^ Jsw 



sw 



1, 



, pQED Jm 



QED 



(60) 
(61) 



With this definition, tan (5 is gauge independent and not affected by QED-corrections 
and infrared divergences. Moreover, the one-loop correction to the decay A^ — > r+r" 
using the renormalization scheme of [1] is quite small [15]. Hence, if tan (5 is defined 
by eq. (l60l . it does not suffer from the numerical instabilities found for the cases of the 
gauge-independent schemes in sections |5] El 



8 Conclusions 

In this paper, the renormalization of tan (3 has been studied in view of the criteria gauge 
independence, process independence, and numerical stability. The DR scheme as well 
as the schemes presented in [1,2] have been shown to imply a gauge dependence of 
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tan P already at the one-loop level. Therefore, three gauge-independent schemes have 
been developed — however, using these schemes produces unacceptably large numeri- 
cal instabilities in higher-order calculations. The conflict between gauge independence 
and numerical stability has been made more explicit by a general statement about a 
large class of process-independent schemes, where tan /3 is defined via quantities of 
the Higgs sector. We have shown that all gauge-independent schemes of this class lead 
to numerical instabilities. Hence, it seems to be impossible to find any renormalization 
prescription that satisfies all three above criteria. 

As a way out of these problems, process-dependent schemes can be used, but such 
schemes also have specific drawbacks. Conceptually, the flavor dependence of tan (3 
is unsatisfactory, and technically, the necessity to calculate three-point functions and 
the possible appearance of QED and QCD corrections are disadvantages. 

In the course of our analysis, two schemes emerge as the best compromises both 
in conceptual and practical respects. On the one hand, the DR scheme is the most ad- 
vantageous among the process-independent schemes. It is technically very convenient, 
numerically perfectly well behaved — and although it is in general gauge dependent, it 
is not gauge dependent at the one-loop level in the context of the important class of 
gauges. On the other hand, defining tan 13 via the decay t~^t~ provides a par- 

ticularly attractive process-dependent alternative. In this scheme, tan (5 is directly con- 
nected to an observable and therefore gauge independent as well as renormalization- 
scale independent. Furthermore, this specific process is theoretically very clean since 
it involves no QCD corrections at the one-loop level and the QED corrections can be 
spht off. 

Both the DR scheme and the t'^t~) scheme have specific advantages, so 

depending on the situation one or the other can be more useful in practice. However, 
finally a decision should be made for one definition of tan /?, since a common renor- 
mahzation scheme is important to allow direct comparisons between different higher- 
order calculations. In the {A^ r+r~) scheme, the specific process is chosen merely 
for technical reasons. From an experimental point of view, the decay A^ t'^t~ is 
only one possibility amongst a variety of potential observables for the determination 
of tan/? (see e.g. [16]); a key observable for the definition of tan/3 does not exist. 
This reflects the fact that tan j3 is an auxiliary parameter. Accordingly, the advantages 
of the process-dependent scheme appear less significant. Owing to its technical con- 
venience and its process independence, we assess the DR scheme as the best choice 
for defining tan j3. 

Acknowledgments: We are grateful to J. Guasch and W. HoUik for valuable discus- 
sions. 

A Slavnov-Taylor operator 

The (gauge part of the) Slavnov-Taylor operator of the MSSM can be written in the 
form 




V='I>1,2,G0 ,aO ,H± , 

other fields 




sr ST 

SY^ Sip 
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^ , sr sr 5r sr y 8t 



(62) 



where "other fields" stands for the (s)quark, (s)lepton, chargino and neutralino fields. 
For the full form also including supersymmetry and translational ghosts, see ref. [4]. 
The fields are auxiliary fields that couple to the gauge-fixing term and are defined 
as the BRS transformation of the Faddeev-Popov antighosts c^. The fields are 
sources for the BRS transformations of the fields (fi. For our purposes, the Feynman 
rules involving ^z,GO,Ao,</)i 2 important. They can be derived from the Lagrangian 
(s denotes the generator of BRS transformations): 

Ly = + Y^^s^2 + YgosG"^ + YaosA^ + Yz^sZ^" + ... 

= y^,\ [9 {ic+{cpG- - spH-) + h.c.) + ~gc^{-cpG'> + s/sA^)] 

+ [9 {-ic-^ispG+ + Cf,H+) + h.c.) + -gc^i-spG'' - 0^^°)] 

+ Yg°\ [9 (c+(-G-) + h.c) + gc^{cp4>i + sp(l>2) + 2c^Mz] 

+ yA^\ [9 {C^i-H-) + h.c.) + 5c^(-s/30i + Cf3^2)] 

+ Yzi^ [di'c^ - gcw{-iW+c- + iWc^)] + ... (63) 



Mm 
Mz • 



with g = \/ g"^ + g'"^ and cw 

B Gauge-fixing terms 

For most parts of the paper we leave the form of the gauge fixing open. However, in 
the one-loop calculations two specific gauge-fixing terms are used. The gauge is 
defined by the choice 

= d^A^, (64) 

J^^ = d^Z^ + Mz^G^, (65) 

:r± = S'^W^ ± iMw^G^, (66) 

Lfix= E 5^t^^ + ||i?^P (67) 

V=A,Z,W+,W- 

for the gauge-fixing functions and the gauge-fixing Lagrangian. For loop calcula- 
tions, it is useful to eliminate the B fields via their equations of motion, yielding 
Lfix = The x-dependent terms (as well as the ghost terms) in the 

Lagrangian that are necessary to satisfy 5(r(°)) = are obtained from Lfix,ghost,x — 
[s + x9^)(c^(jF^ + where s denotes the BRS operator. For ehminated B 

fields they read 



-gauge _ 1 



c'^d^A^ + c^id^'Z^ - ^MzG°) 



+ c-{d^'W+ - i^MwG+) + c + idf^W- + i^MwG') 



(68) 
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In the generalized gauge defined by eq. d30l . the correct form of the x"^"^" -dependent 
terms is given by 

L^zAO = x^^^c^MzA^ (69) 

C General gauge-independent scheme 

In this section we show that the most general gauge-independent renormalization 
scheme of the class ( 15^ is given by eq. ( I5^ . 

Since the 771,3 scheme is gauge independent, any other gauge-independent scheme 
has to satisfy 

di,{5tp-{5tp)m,) = Q, (70) 

leading us to the question on which combinations of the quantities in eq. ( 15^ are 
gauge independent. In order to answer this question it is sufficient to study linear com- 
binations of the quantities in eq. (l52l . leaving away S^o^o f^f^ and S^o^ since the 
gauge-independent combinations Ka^h.k involving S^o^o are already known 

and since S^o^ itself is a linear combination of Tij^oqo and 5ti 2 as a consequence of 
the Slavnov-Taylor identity. 

The gauge dependence of the remaining quantities can be expressed using the iden- 
tity 5(r(i)'^^s) = as 

{6ti,6t2, T.ji^aQO,T,Hh{MH), T,Hh{Ml), S^o^o, S'^//, ^'hhV 

- A/f /'r(l)''^<=S T.(l),reg p(l),reg p(0) y.(l),reg ^(0) 

r(l).i'eg/,^2 ^, r^^)''^^S( ]\/f2\ -p(l),reg p(l),reg -p(l),reg\'^ 

i XHY^ K^^-^h),^ XhYH ^^^-^h), >- ^AOY^o ' ^ xHY„^^ xhY^ j ' ^'^) 

where M denotes an invertible matrix whose entries consist of tree-level expressions 
likeri°|, etc. 

Calculation of the vertex functions involving x shows that in the gauge there is 
only one vanishing linear combination, namely spT^^y^'^^ ~ '^P'^xY^'^^ ~ ^' "^^^ other 
vertex functions involving x in e<l- (ItTI are linearly independent. Correspondingly, the 
combination d29l is the only gauge-independent linear combination of the quantities 
on the l.h.s. of eq. (fTTT i in the i?^ gauge. However, in more general gauges this combi- 
nation is not gauge independent (see eq. (l34l ). Hence, there is no gauge-independent 
linear combination of the quantities on the l.h.s. of eq. (f/TT i. 

This confirms that the expressions KA,H,h are the only gauge-independent combi- 
nations of the quantities in eq. i57i . as was to be shown. 
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